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ABSTRACT 

We present the kinematic results from our ARGOS spectroscopic survey of the Galactic 
bulge of the Milky Way. Our aim is to understand the formation of the Galactic bulge. We ex- 
amine the kinematics of about 17,400 stars in the bulge located within 3.5 kpc of the Galactic 
centre, identified from the 28,000-star ARGOS survey. We aim to determine if the formation 
of the bulge has been internally driven from disk instabiUties as suggested by its boxy shape, 
or if mergers have played a significant role as expected from ACDM simulations. From our 
velocity measurements across latitudes b = -5°, -7.5° and -10° we find the bulge to be a 
cylindrically rotating system that transitions smoothly out into the disk. Within the bulge, we 
find a kinematically distinct metal -poor population ([Fe/H ] < —1.0) that is not rotating cylin- 
drically. The 5% of our stars with [Fe/H] < —1.0 are a slowly rotating spheroidal population, 
which we believe are stars of the metal-weak thick disk and halo which presently lie in the in- 
ner Galaxy. The kinematics of the two bulge components that we identified in ARGOS paper 
III (mean [Fe/H] « -0.25 and [Fe/H] w H-0.15, respectively) demonstrate that they are likely 
to share a common formation origin and are distinct from the more metal-poor populations of 
the thick disk and halo which are co-located inside the bulge. We do not exclude an underlying 
merger generated bulge component but our results favour bulge formation from instabiUties 
in the early thin disk. 

Key words: Milky Way Galaxy - Galaxies - Stellar Populations - Kinematics and Dynamics 
- Spectroscopy - Galactic Bulge. 



1 INTRODUCTION 

The principal goal of the ARGOS spectroscopic survey is to obtain 
chemical and kinematic data for a substantial number of stars to un- 
derstand the origin of the Galactic bulge of the Milky Way. The two 
main scenarios for bulge formation are hierarchical merging as seen 
in Lam bda Cold Dark Ma tter (ACDM) simula tions of Galaxy for- 
mation i Abadi et al .l2003h and disk nstabi lity jCombes & SandersI 
ll98lURaha et al.ll99ll ; lAthanassoulal2005l) . These different forma- 
tion processes imprint distinct kinematic signatures and thus can be 
tested using observations of the stellar kinematics of the bulge. 

Boxy/peanut bulges rotate almost as rapidly at high latitude 
as in the plane (cylindrical rotation). This has been seen in obser- 
vations of boxy bulges in other galaxies, e.g. lFalcon-Barroso et al.l 
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( l2006l) , and in numerical simulations ( I Athanassoulall20o"2h . This is 
unlike the merger-gener ated 'classical' bulges which show slow ro- 
tation at high latitudes. [Howard et al.l ( |2009|) have identified such 
cylindrical rotation in the Galactic bulge for a subsample of the 
more metal-rich stars. Shen et al. (2010) constrain the bulge to have 
no underlying merger-generated component to a level of ^5 8%. 

The ARGOS survey sample is unbiased in metallicity, and 
provides kinematic data of stars in the Galactic bulge over the 
wide range of abundance that is a consequence of the bulge for- 
mation process. This may enable us to investigate at what stage the 
boxy bulge structure was formed during the evolution of the inner 
Galaxy. 

We test the prediction of rapid high-latitude rotation by mea- 
suring velocities for bulge stars in 25 fields at three Galactic lati- 
tudes in the south (plus three fields in the northern bulge as checks 
on symmetry), out into the disk across a longitude range from +26° 
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to —31°. Our sample size is large enough to detect a 5% contribu- 
tion from a classical slowly-rotating merger-generated spheroidal 
component un derlying the dominant boxy/peanut bu l ge. We note, 
however, that ISaha. Martinez- Valpuesta. & GerhardI (1201 2|) have 
shown that a weak classical bulge would be spun up into rapid 
rotation by the boxy/peanut bulge, and would therefore be diffi- 
cult to detect kinematically. The ARGOS survey would also in- 
clude stars of the metal-poor "first stars" population, if they are 
present in the inner bulge as predic ted by cosmological simulations 
(iDiemand. Madau. & Moorell 20051) . 

This paper first reintroduces, in Section 2, the s tellar pop- 
ulatio ns of the bulge identified in ARGOS paper III jNess et al] 
l2013h . Sections 3-5 briefly recapitulate the details of the obser- 
vations, stellar sample selection and the methods used to measure 
stellar parameters and determine distances. The distance determi- 
nation has enabled selection of a relatively uncontaminated sample 
of stars in the inner Galaxy. Using our estimates of stellar param- 
eters, we have been able to identify and eliminate foreground and 
background stars. In Section 6, we present our results on the cylin- 
drical rotation profile obtained for the bulge and the slow rotation 
of the metal-poor population with [Fe/H] < -1.0 in the inner re- 
gion of the Galaxy. We also compare the kinematics of populations 
of different metallicity and provide evidence for a common dynam- 
ical history for stars with [Fe/H] > —0.5. In Section 7, we compare 
our kinematics with those for N-body models and discuss the origin 
of the components we find in the bulge in Section 8. Conclusions 
are presented in Section 9. 



2 STELLAR POPULATIONS IN THE BULGE 

ARGOS paper 111 jNess etal]|2013h presents the metallicity distri- 
bution of the bulge. We find that the bulge is a composite popula- 
tion. Figure [T]shows our decomposition of the ARGOS metallicity 
distribution functions (MDFs) at 6 = —5°, —7.5° and —10° for 
stars in the bulge with Rq < 3.5 kpc across I = ±15°, into five 
Gaussian components. These components are labeled A-E in or- 
der of decreasing [Fe/H]. We used Gaussians for simplicity in this 
analysis and the number of components fit was motivated by the 
distinct peaks seen in the metallicity distribution function and a 
statistical analysis; components A-C are our primary components 
and D and E are included to represent the small fraction of very 
metal poor stars in the sample wh ich increase in fr action further 
away from the Galactic plane (see iNess et al.ll2013l for a discus- 
sion). These components are proxies for stellar populations and are 
used as a tool to measure the changing contribution fraction of stars 
with different metallicities as a function of (I, b). The analysis of 
the m etallicity-depende nt split red clump stars along the minor axis 
fields jNess et alj|2012h has demonstrated that stars with [Fe/H] > 
-0.5 are part of the boxy/peanut bulge and stars with [Fe/H] < -0.5 
are not. 

If the formation of the bulge of our Galaxy i s indeed a re- 
sult of a dynamical in stability of the initial disk jAthanassoulal 
l2005l : ISheneTal]|2010h then it is critical to interpret and anal- 
yse the bulge in terms of the contribution of these stellar com- 
ponents which reflect the way in which the bulge was fomed 
from the pre-existing disk. For the bulge stars with [Fe/H] > -0.5, 
there are two populations (components A and B). Component A 
is a metal-rich, relatively thin boxy/peanut-bulge population with 
a mean [Fe/H] = -nO. 15 and is concentrated close to the plane. We 
associate this component chemically with the thin disk in the in- 
ner region. Component B is a more metal-poor (mean [Fe/H] = 
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Figure 1. MDFs for stars within Rq < 3.5 kpc at from left to right b = 
-5°, b = -7.5° and b = -10°, for I = ±15°, showing the changing 
contribution of metallicity fractions with latitude. 

-0.25) and thicker boxy/peanut-bulge structure, distributed across 
6 = -5° to 6 = -10°. We have argued in ARGOS 111 that com- 
ponent B represents the stars formed out of the thin disk at early 
times. 

In summary, contributing to the stellar density in the bulge 
region are the thin boxy/peanut-bulge (A), the thick boxy/peanut- 
bulge (B), the thick disk (C), metal-weak thick disk (D) and halo 
(E), in order of decreasing [Fe/H]. The stars of the halo are probably 
not part of the bulge, in that they likely have a different origin from 
the majority of the stars in the bulge. These halo stars are likely on 
highly eccentric orbits passing through this region at this time. 



3 OBSERVATIONS AND STRATEGY 

We acquired our spectra with the fiber-fed AAOmega system on the 
Anglo Australian Telescope. Figure|2]shows the locations of the 28 
fields, each of two-degree diameter, selected for this program. We 
observed about 28,000 stars between 2008 and 2011. Our survey 
covers the bulge, which is estimated from the 2MASS stellar profile 
to cov er a longitude rang e of ±15° and extend in latitude to about 
±12° jPwek et alJI 19951) . Our fields extend in longitude out into 
the thin and thick disk, in order to measure the stellar kinematics at 
the transition of the bulge and disks. 

Each field contains about 1000 stars and a full observation at 
each {I, b) requires 3 separate setups of the AAOmega fibres. About 
350 stars can be observed simultaneously, and each field of 1000 
stars was split into three magnitude intervals, as a partial proxy 
for distance along the line of sight through the bulge. In our setup 
we used 25 sky fibres for sky subtraction, and 8 fiducial fibres for 
tracking alignment. 

The stellar parameters were derived from a comparison 
with synthetic spectra generated in the ID Lo cal Thermal E quilib- 
rium (LTE) stellar synthesis program M OOG (|Snedenlll973l) usin g 
the Castelli/Kurucz model atmospheres jCastelli & Kuruczll2004l) . 
From the distance estimates for our stars, which are based on the 
stellar parameters, we have isolated the 17,400 stars in our sample 
within y — ±3.5 kpc (here j/ is a Cartesian coordinate along the 
centre-Sun line, with origin at the Galactic centre). 

The full details of our surv ey, the observations a nd the anal- 
ysis, are provided in ARGOS II jFreeman et al ] |2013h . The metal- 
licity decomposi tion of populations are discussed in ARGOS III 
jNess et al.l20T3h . 



4 MEASURING RADIAL VELOCITIES 

The radial velocities of stars in our fields were measured from our 
medium resolution spectra (R = 11, 000) covering the spectral re- 
gion 8000 — 88OOA which includes the strong Ca II infrared triplet 
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Figure 2. Galactic latitude and longitude for the 28 two-degree-diameter 
fields in our survey. The filled circles indicate fields that were completely 
observed. We are missing stars only in one field, at {I, b) = (20°, —7.5°) 
for which only 600 star's were observed (rather than 1000 typically per 
field). The itmer rectangle represents the bulge region. 



absorption lines. Following the data reduction and sky subtraction, 
the data were cleaned using a sigma clipping algorithm to remove 
residual sky lines which improved the accuracy of radial velocity 
measurements. Radial velocities were obtained via cross correla- 
tion of each observed spectrum with a number of template stars 
generated from synthetic spectral models. 

These stellar parameters of the template stars covered a wide 
range of [Fe/H] at a fixed gravity and temperature. The metallicity 
values used were [Fe/H] = -2.0, -1.0 and 0.0, for Te// = 5000K 
and log 17= 2.8 (as we ex pect most of our stars to be clump giants 
& Maol200ll) ). The best matching template was chosen 
for each star to measure the radial velocity, as determined by the 
smallest error of the fit in the cross correlation. The error at this 
resolution and signal to noise (about 50 to 80) is f» 0.9 km s"^ (see 
ARGOS II). All heliocentric radial velocities were transformed to 
Galactocentric velocities (see Section 5). 



5 BULGE MEMBERSHIP 

From the ARGOS spectra, we have measured radial velocities, stel- 
lar parameters ( Teff, log (7) and chemical abundance data ([Fe/H], 
[a/Fe]) for our stars. The stellar parameters allow us to identify 
foreground dwarf stars and, from the estimated distance, to deter- 
mine which stars are foreground and background giants. 

The rotation and dispersion results (see Section 6) are for stars 
chosen to lie within y = ±3.5 kpc where y is a Cartesian coordi- 
nate along the centre-Sun line, with origin at the Galactic centre). 

This selection on y is made, rather than a cut in Galactocen- 
tric radius Rq ^ 3.5 kpc, so as to include the stars in our fields 
at Galactic longitudes |/| > 26°, which have a minimum galacto- 
centric radius, Rq > 3.5 kpc. In total 17,400 stars remain after 
this distance cut, about 70% percent of the original sample. The re- 
maining 30% of the stars are in the foreground or background. For 
studying the properties of the bulge, it is important to remove these 
contaminants. 



6 RESULTS 

6.1 Cylindrical rotation 

Our rotation curves are given in the left panel of Figure|3] and show 
the Galactocentric rotation velocity in all of our fields out to longi- 
tudes of —30° to -|-26°, made using a cut of \y\ < 3.5 kpc, for stars 
with [Fe/H] > —1.0. We find cylindrical rotation for the bulge, 
with only a weak dependence of rotation speed on latitude. These 
data extend the cylindrical rota tion fou nd in the BRAVA survey 
(Howard et al. 2009; Kunder et alllOlj) to longitudes of |;| = 10°. 
The BRAVA data jKunder et alj20I2r) are included for comparison 
in Figure[3]at their latitudes of & = —4°, —6° and —8°. 

We corrected for the solar reflex motion, adopting the 
local standard o f rest velocity at the Sun to be 220 km s~^ 
jKerr & Lvnden-B ell 1986) and a solar peculi ar velocity of 16.5 
km s ~^ in the direction {l,b) = (53°, 25°) jMihalas & BinnevI 
Il98lh . The corrected velocity, the Galactocentric velocity Vgc is 
then 

Vgc = Vhc + 220[sin(0 cos(fe)] -f 16.5[sin(6) sin(25) + 
cos(6) cos(25) cos{l — 53)] 

where Vhc is the heliocentric radial velocity in km s^^ and angles 
(l,b) are in degrees. 

Each data point in Figure |3] represents the mean velocity of 
about 600 ARGOS stars. The error bars show the sampling error in 
each field, but due to the large size of our sample they are, in most 
cases, smaller than the symbol size. Note that the mean velocities 
of our fields along the minor axis is non-zero: the mean velocities 
for our three minor axis fields are between about -6 km s^^ to -1 1 
km s^^ . These non-zero values may come from streaming motions 
in the bulge, or our adoption of a Sun-centre distance of 8 kpc. 
Alternatively the offsets could be a sampling effect derived from 
preferential sampling of more metal-rich stars on the near side of 
the bulge, which are closer to the plane. Section 6.3 demonstrates 
that the metal rich stars have lower mean velocities in our fields 
than the metal poor stars across our longitude range. The rotation is 
faster on the far side of the bulge at negative longitudes. This may 
come about because, on the far side of the bulge, we are observing 
further along the bulge/bar at a given given longitude, relative to 
the near side. Systematic distance errors may also contribute. The 
rotation increases to > 150 km s~^ at our largest longitudes, which 
are in the Galactic disk at / = —31° and +26°. 

The results for our latitudes of 6 = -5°, -7.5°, b = -10° 
and also a few fields at —4° and +10° are tabulated in Table [T] 
Although we do not share any common fields, our velocity results 
compare fairly well with those from the BRAVA survey o f M giants 
in the bulge (see lHoward et al.l2009l : lKunder et alj2012l') . from 6 = 
—4° to 6 = —8° and mostly across I = ±10°. As the lowest 
latitudes of the BRAVA fields are closer to the Galactic plane, we 
expect (and confirm) a slightly faster rotation for the stars in their 
b = —4° field compared to our lowest latitude fields b — —5°. 
The velocities at positive longitudes at 6 = —4° from BRAVA are 
higher than the ARGOS measurements at b = — 5° which are better 
matched by the = —6° BRAVA fields. The rotation and dispersion 
shows symmetry about the major axis, as seen by comparing the 
fields at negative and positive latitudes. 

The rotation curve shown in Figure |3] is simply the average 
value of Vgc as a function of longitude, uncorrected for projection. 
In this representation, of rotation as a function of longitude, the 
true rotation curve for the Galaxy would be the mean value for 
the azimuthal component V^ of the stellar velocity for stars with 
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Figure 3. Rotation curves (left) and velocity dispersions (right) for the 16,600 stars in our survey in the bulge region at y = it 3.5 kpc from the Galactic centre 
with [Fe/H] > -1.0. The red diamonds are 6 = —5°, yellow stars are h = —7.5°, blue rectangles are b = —10° and black squares are h = +10° ARGOS 
fields. The two curves correspond to our data at b = —5° and b = —10 °. The horizontal error bars represent the two-degree size of the fields and the vertical 
error bars are the sampling errors. The BRAVA data dKunderet al.l2012h are also shown using the smaller syinbols which correspond to their fields at 6 = —4° 
,6 = -6° and 6 = 8°. 
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Figure 4. The two versions of the rotation curve for the model at fe = —5° 
as a function of longitude. The red points are calculated as the mean Vqc 
as for the data in Figure |3] The black points show the azimuthally aver- 
aged V0 component of the stellar velocity at the minimum value of Rq at 
each longitude, i.e. the true rotation. The red points from the model are for 
comparison to the data at 6 = —5° in Figure[3] 



bulge/bar lies at an angle of about 20° to the Sun-center line, 
fields at negative longitudes are further from the Galactic centre 
than fields at the same |Z| on the positive side. All three ranges of 
magnitude are included in each data point, except for our field at 
{I, b) = (20°, —7.5°), which is missing the faintest stars. The ve- 
locity dispersion in the bulge is higher at lower latitude, most sig- 
nificantly in the central region, and the dependence of dispersion 
on latitude decreases with Across longitude, the dispersion at 
& = -10° is flatter than at & = -5°. 

The dispersion plot in the right hand panel of Figure|3]includes 
the data from the B RAVA survey at 6 = -4°, -6° and h = -8° 
jKunder et al]|2012l) . The dispersion measured from the BRAVA 
survey is similar to the ARGOS results. The h — —8° fields from 
BRAVA have dispersions at about the same level as the —10° AR- 
GOS fields rather than the higher dispersion as measured in the 
ARGOS b = -7.5° fields. In the inner region |i| < 5°, the disper- 
sion is about 20% higher in the BRAVA fields at —4° compared to 
the ARGOS fields at —5°. As discussed in paper III, there appears 
to be a structural change in the bulge populations in the bulge at 
b < —5° and it appears that this is most significant within [/| < 5°. 
The velocity dispersion is metallicity sensitive and is lower for the 
more metal-rich populations. This is examined in Section 6.3. 



an Rq value close to the minimum Rq at a given longitude. We 
can use our N-body model as a guide to the difference. Figure |4] 
sh ows the two rotation curves as derived for the model that we used 
in iNess et alj (|2012!). The red points show the mean Vbc which 
corresponds to the curve shown in Figure|3]at a latitude of & = — 5° . 
The black points in Figure |4] show the mean as a function of 
longitude. At \l\ > 10° the difference between the mean Vbc and 
the mean is only a few percent, but in the inner regions the V,/, 
curve is steeper. The red points in Figure |4]represent the data well 
at 6 = —5° in Figure|3] 

The line-of-sight velocity dispersion across our fields is shown 
in the right hand panel of Figure |3]for stars with [Fe/H] > —1.0. 
Note that no a-clipping has been done for these measurements 
which are the standard deviation of all stars with [Fe/H] > —1.0 
about their mean velocity in each field. As the major axis of the 



6.2 Rotation of the metal-poor stars 

We now proceed to investigate the rotation of the sample of stars 
with [Fe/H] ^ —1.0. These stars comprise about 5% of the data, or 
about 750 bulge stars within a distance \y\ < 3.5 kpc. The rotation 
curve and dispersion profile for all stars with \y\ < 3.5 kpc and 
[Fe/H] < —1.0 are shown in the left and right panel of Figure |5] 
respectively. 

The mean rotation for the metal-poor stars is lower than for 
the metal-rich stars ([Fe/H] > —1.0), by about 50%. The rotation 
in both high (b = —10°) and low (b — —5°) latitude fields is sim- 
ilar at positive longitudes on the near side of the bulge. On the far 
side of the bulge the high latitude dispersion has a flatter profile. 
Note that the sample size in some fields is small; there are only two 
stars in the field at {I, b) = ( — 20° , —5°). The low metallicity stars 
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Figure 5. Rotation (left) and dispersion measurements (right) for tlie 750 stars within \y\ < 3.5 kpc with [Fe/H] -1.0. The red diamonds are for b = —5°, 
the yellow stars are 6 = —7.5°, the blue rectangles are 6 = —10° and the black circles are for b = +10°. The error bars in the vertical direction represent 
the sampling enor for each field. 



show a net rotation, and this rotation persists to the lowest metal- 
licities [Fe/H] < —1.5, but the number of stars is small and the 
error bars correspondingly larger (see Figure 0. There is a larger 
fraction of more metal-poor stars at high latitudes than low lati- 
tudes, with about 400 stars al b — —10°, 100 stars in the three 
fields at b — +10° and 100 stars in ou r fields at 6 = —5° and 
6 = -7.5°. lHarding & Morrisonl ( Il993l) had earlier noted a simi- 
lar change of kinematics with metallicity in a single bulge field at 
(/,6) = (-10°,-10°). 

The velocity dispersion of the metal-poor stars in Figure |5]is 
significantly higher than that of the metal-rich stars in Figure [3] 
Metal-poor stars in fields at a latitude of 6 = —10° (blue symbols) 
show a greater increase in dispersion than the fields at 6 = —5° 
(red symbols), relative to stars at [Fe/H] > -1.0. At fe = ±10°, 
the dispersion on the far side of the bulge is higher than the near 
side of the bulge, shown in the left panel of Figure |5] No sigma 
clipping has been used in the plot, because the sample is small and 
we did not want to exclude high velocity outliers. No high velocity 
outliers are seen at this latitude at positive longitudes. At negative 
longitudes there is a group of stars at 6 = —10° with velocities 
> 300 km s~^(see the right hand panel of Figure|8). Four of the 3 1 
stars in the field (/, b) = (-20°, -10°), with an [Fe/H] of -1.2, 
have velocities between 300 — 330 km s^^and lie within Rq < 
4.0 kpc of the Galactic centre. If these four stars are removed, the 
velocity dispersion decreases from a = 175 km s~^ to a = 125 
km s~^. The higher dispersion seen for more metal-poor stars is 
discussed further in Section 7.1. The asymmetry in the dispersion 
seen for the stars at 6 = —10° in Figure |5] also suggests there 
may be a relationship between the halo stars of the inner Galaxy 
and the bar. In simulations which include a stellar halo component, 
the halo in the inner part becomes elongated and makes a so called 
halobar. Assuming that a non-negligible number of stars follows the 
dynamics of the inner part of the halo, this would provide a metal- 
poor compo nent of halo stars tha t may resemble the dynamics in 
Figurel5l(see lAthanasso"uial ( |2007|) ). 

6.3 Kinematics of stars with [Fe/H] > 1.0 

In paper II we identified stars with [Fe/H] > -0.5 to be of the 
bulge, and we believe these stars had a common dynamical his- 
tory. We decomposed our sample of observed stars into metallicity 
components representative of the populations in the inner region. 



The most metal-rich component A with mean [Fe/H] ~ +0.15 is a 
relatively thin and centrally concentrated part of the boxy/peanut- 
bulge. Component B with mean [Fe/H] ~ —0.25 is a thicker 
boxy/peanut-bulge structure seen in relatively constant fraction 
across our latitude range in the inner region. We speculate that the 
stars in component C, with mean [Fe/H] ~ —0.70, are from the 
early thick disk which may have been puffed up in the inner region 
during a bar buckling event. We note again that component C does 
not appear to be part of the boxy/peanut structure of the bulge. The 
components are overlapping in metallicity, and component B ex- 
tends down to near [Fe/H] = —1.0, well into the metallicity range 
of component C. Therefore, stars as metal-poor as [Fe/H] ~ —1.0 
are all contributing to the bulge in the inner region. The stars with 
[Fe/H] < — 1.0 show slow rotation and high dispersion (see Figure 
O and we attribute them to the metal-weak thick disk and halo. 

We examine the kinematics of stars in metallicity intervals 
coiTesponding to our components A, B and C to demonstrate that 
A and B share common properties distinct from C. Figure |6]shows 
the rotation (top panels) and velocity dispersion (bottom panels) 
for stars in metallicity intervals [Fe/H] > (A), —0.5 < [Fe/H] 
s; (B) and -1.0 < [Fe/H] sC -0.5 (C). The rotation for all 
three [Fe/H] intervals is fairly similar, with component B showing 
the fastest rotation, about 20% faster than components A and C at 
b — —5°. Although the rotation curves are similar, the velocity dis- 
persion profiles are clearly different between the three components. 
The dispersion profiles for components A and B have similar mor- 
phologies with changing latitude; component A is a kinematically 
colder replica of component B. We note that component A shows a 
relatively high dispersion at 6 = ±10° and I = —5°, +10° which is 
seen to a lesser extent also in component B. These fields correspond 
to similar distances along the bar of about 2.0 kpc. The dispersion 
for the most metal-rich fraction of stars (A) is about 30% lower 
than for stars associated with component B. The relative change of 
dispersion with latitude for components A and B is consistent with 
a lower scale height for component A than for B. 

Component C however has a clearly different dispersion pro- 
file and latitude dependence relative to stars with [Fe/H] > —0.5, 
with a relatively latitude-independent dispersion out to |Z| = 10°. 
At larger longitudes, |Z| > 10°, the dispersion decreases out into 
the thick disk further from the Galactic centre (see the far right 
panel in Figure|6l. The common rotation shape and dispersion pro- 
files seen for components A and B suggest that A and B share a 
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Figure 6. Rotation (top panel) and velocity dispersion (bottom panel) across the bulge to |/| = 20° for the 16,600 stars with [Fe/H] > —1.0 within < 3.5kpc 
of the Galactic centre. The three plots correspond to different metallicity bins, from left to right in decreasing [Fe/H] as shown. Note that the discrete bins 
are used to represent stars of our components A,B and C from left to right. Although the rotation curves are similar, the dispersion clearly demonstrates the 
difference in kinematics of stars with [Fe/H] > —0.5 and with [Fe/H] < —0.5. There are 3100, 8600 and 4900 stars in each plot, from left to right. The red 
diamonds are for b = —5°, the yellow stars are b = —7.5°, the blue rectangles are b = —10° and the black circles are for b = +10°. 



similar formation history or common formation mechanism. The 
similar rotation but different dispersion for component C indicates 
that it is a distinct population, although it is part of the central ro- 
tating bulge in the inner region. This could be due to the different 
initial velocity dispersion of the thin and thick disks before the in- 
stability event. 

6.4 Stellar streams in the bulge 

The top panel of Figure |7] shows the rotational profile of the bulge 
as a function of [Fe/H] in the five bins of positive longitude, where 
the stars are divided into metallicity bins of 0.5 dex, except for the 
most metal-poor stars of which there are very few, so the lowest 
[Fe/H] bin takes all stars from the most metal-poor star identified 
([Fe/H] = —2.8) to [Fe/H] = — 1.5. The most metal-rich bin is from 
[Fe/H] = to 0.5 dex. These bins correspond approximately to the 
components which we label A-E, in order of decreasing [Fe/H], as 
outlined in Section 2. 

The upper panel of Figure |7] shows the existence of a group 
of metal-poor counter-rotating stars in the bulge at positive longi- 
tudes. The negative longitudes examined in this way show the same 
overall kinematic behaviour as a function of [Fe/H] with latitude, 
but do not reveal any distinct counter rotating-groups in any [Fe/H] 



bin. The counter-rotating stars are shown in the shaded panels in 
the fields (/,6) = (10°, -7.5°) and (/, 6) = (10°, -10°) . There 
are seven stars in total in these fields with [Fe/H] < — 1.5 and these 
are shown individually in Figure |8] Six of these stars which are 
shown in the left hand panel of Figure [8] are also co-located in dis- 
tance at around 4.9 kpc along the line of sight from the Sun, with 
a distance spread of a = 0.6 kpc (see Figure |9}. These six stars 
have a metallicity of (/i([Fe/H]), cr([Fe/H])) = (-1.68, 0.03) and a- 
enhancement range of (/i([a/Fe]), (7([a/Fe])) = (0.50,0.17) and the 
five stars in the field at {I, b) — (10°, —7.5°) are moving with ve- 
locities of between —90 and —133 km s~^. There is one star in 
this group located in the field {I, b) = (10°, -10°), with [Fe/H] = 
-1.6andl/GC= -148 km s"\ This group may be an independent 
stellar stream in the bulge at a distance of around 4.9 kpc from the 
Sun, with a line of sight velocity of (^(Vgc), cCVgc)) = (-115 km 
s~^, 22 km s~^). From its location it seems this group could be 
related to the comers of the bar. Comparing the percentage of stars 
in this group by using their K-magnitude compared to the 2IV1ASS 
stars in the field, about 200 stars are expected to belong to this 
low-metallicity group from of all stars in the 2MASS count for this 
field. 

The right hand panel of Figure|8]shows four stars at negative 
longitudes which are tightly clumped in distance (as noted in Sec- 
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Table 1. Galactocentric velocities observed in our fields across the bulge 
for stars with [Fe/H] > -1.0. 



tion 6.2) and velocity, with (^i(VGc), (^(Vgc)) = (320 km 8 
km s~^). This group of stars lies in the disk and is located at about 
8.4 kpc away from the Sun along the line of sight with a distance 
spread of a = 0.34 kpc (see Figure |9}. These stars are also very 
similar in their metallicity, with (^([Fe/H]), a([Fe/H])) = (-1.24, 
0.09) and have an a-enhancement range of (^([a/Fe]), a([a/Fe:])) 
= (0.38,0.15). Comparing the percentage of stars in this group by 
using their K-magnitude compared to the 2MASS stars, a lower 
limit of about 100 stars are expected to belong to this group from 
all stars in the 2MASS count for this field (not taking into account 
2MASS incompleteness given their mean K-magnitude of about K 
= 13.2). 

The panels in Figure |7] span longitudes to the edge of the 
bulge, at I = 0°, I = 5° and I — 10° and also our fields out- 
side the bulge out into the disk at I — 15° and I — 20°. Three 
latitudes zones are included in these plots, at 6 = —5°, —7.5° and 
— 10° . Missing points in the plots indicate no stars in this [Fe/H] bin 
for the corresponding (/, 6). From these rotation plots as a function 



of [Fe/H], it is clear that the thicker boxy/peanut bulge component 
B, corresponding to > [Fe/H] > -0.5 is generally the fastest ro- 
tating. This component rotates most rapidly in the lower latitude 
fields, b = —5°. The rotation speed increases with longitude. The 
increase is more rapid for the metal-rich components, at around 48 
km s^^kpc^^ for [Fe/H] > —1.0, compared to about 27 km s~^ 
kpc^^ for the metal-poor stars with [Fe/H] < —1.0. 

All stars with [Fe/H] > —1.0 (components A,B,C) have fairly 
consistent relative rotation trends across longitude, and components 
A and B show the largest difference in rotation velocity with lati- 
tude near the comers of the bulge, at I — 15°. The trends seen 
with longitude in the more metal-poor stars ([Fe/H] < —1.0) are 
not as consistent as for the more metal-rich ([Fe/H]> —1.0) stars. 
Component D rotates more slowly than A,B,C; the rotation of com- 
ponent E depends on the field. 

From the bottom panel of Figure |7] it is clear that the veloc- 
ity dispersion mostly increases as metallicity decreases down to 
[Fe/H] < —1.0 ; for lower metallicities it can decrease or increase 
depending on the field. Note that in the field where the maximum 
dispersion is found, a ~ 215 km s^^ at {I, b) — (0°, —10°), there 
are only 9 stars in this bin of [Fe/H], one of which has a Vgc = 
412 km s~^ (there is no sigma clipping: all stars are included in 
the plot). The error bars become larger across all longitudes for the 
metal-poor stars, [Fe/H] < —1.0, as the number of metal-poor stars 
is small. 



7 COMPARISONS TO MODELS 

In this section we compare our data to a series of six simulations 
of a rotating disk embedded in a live dark matter halo. None of our 
models contains gas, i.e. there is no star formation and no informa- 
tion on metallicity for us to compare to. We also include only one 
disk component, instead of two separate thin and thick disk com- 
ponents. Th e disk and the halo in itial density distributions are as 
described in lAthanassoulal ( |2007|) and have been built in equilib- 
rium in each other's potential as described in that paper. All models 
have identical initial conditions for the disk, with a Toomre param- 
eter equal to Q=1.2 and a vertical height of 0.2 times the disc scale 
length. In this way some effect from the thick disk is included. We 
consider halo core sizes in the range between 0. 1 and 4 times the 
initial disc scale length. The former leads to model with a submaxi- 
mum disk, while in the latter the disc is maximum (for a discussion 
on the effects of maximum and sub-maximum disks on bar growth 
and evolution, see lAthanassoula & Misirioti3 l l2002h. Their evolu- 
tion w as followed using the Gyrfalcon N-body code ( lDehnenl200Cl 
I2OO2I) . 

All models give rotation curves that compare reasonably well 
with the data, but the changes in velocity dispersion across lon- 
gitude and offsets between b — —5° and b = —10° are model- 
dependent. Although the kinematic profiles of these models are 
similar to those observed, none fits across all latitudes. The vari- 
ation of velocity dispersion across longitude is sensitive to the bar 
properties. For bars that are too short with respect to the observa- 
tions, the dispersion decreases more rapidly than observed at the 
highest latitude and across the longitude range \l\ < 20°. The dis- 
persion profile at the lowest latitude is sensitive to the length of the 
bar which offers a complementary test of the observed scale length 
of the bulge. 

As already discussed in Section 2, we have associated stars 
with [Fe/H] > -0.5 with the boxy/peanut bulge. Our model shows 
the X-shape structure in the density distribution of the stars at high 
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Figure 7. The rotation (top) and dispersion (bottom) versus [Fe/H] in 0.5 dex bins across positive longitudes from Z = 0° to i = 20° . Diamonds are for 
b = —5°, stars are b = —7.5°, rectangles are b = —10°. Shaded panels indicate a group of stars in the inner regions which are at large negative velocities 
and low metallicities. 



latitudes in tlie inner region, showing a very good match between 
observational data for stars with [Fe/H ] > -0.5 at the minor axis, 
similarly to the model discussed in iNess et al ] ( l2012h . Figure 1101 
compares the kinematic data for the stars with [Fe/H] > -0.5 to 
the model. The model velocities have been scaled in amplitude to 
best match the dispersion at 6 = —7.5° out to \l\ = 20°. The rota- 
tion curves of the model (top panel of Figure [TOt are a reasonable 
fit to the data for the stars with [Fe/H] > -0.5. In both the model 
and observations, at I — i: 20° the rotation velocity at fe = —5° 
is approaching 150 km s^^ and at 6 = —10° it is about 120 km 
s^^. The rotation velocity is marginally high in the model com- 
pared to data at positive longitudes at fe = —10°. In the model, 
the Sun is located at 8 kpc from the Galactic centre along the line 
joining the Sun to the centre of the Galaxy and the bar is point- 
ing 25 degrees into the first quadrant. Adjusting the bar angle by 
±10° relative to the Sun-centre line changes the gradient of the ro- 
tation by about 15%: larger bar angles give flatter rotation curves. 
Although the dispersion profile is model dependent, the character- 
istics of the profile across (l, b) are generic to our N-body models 
of instability-generated boxy/bulges. 

From the bottom panel of Figure [TO] it is clear that the over- 
all shape of the dispersion profile in the model matches well the 



observations across our fields for our more metal-rich stars (com- 
ponents A and B, see Figure |6j. At fe = —5°, the observed and 
model dispersions are a fairly good match except at I < —10° 
where the model dispersion decreases faster than the observations. 
At fe = —10°, the observed and model dispersion are both fairly 
flat across the longitude range but in the inner region the obser- 
vations are about 15% lower than the model. The gradient of the 
model and data are, however, similar at the three latitudes, and at 
fe — —5° both model and data show an asymmetry at Z = ± 5°, 
where the dispersion at positive longitudes is higher than the disper- 
sion at the negative longitudes. The feature in the dispersion seen 
near the edges of the bar at I = —5° and I = —10° at the higher 
latitude of fe = —10° (see also Figure |6}, is not reproduced in the 
model. 

In Figure [TT] we compare the model to stars with [Fe/H] > 
-1.0. Although the dispersion of stars in our component C does 
not reflect that of the mode l and is not involved in the boxy/peanut 
structure dNess etaLllioii) . we include this comparison because 
this component is still cylindrically rotating similarly to A and B. It 
is possible that there are a smaller number of stars within the model 
contributing to this dispersion distribution seen in C. We cannot 
break up our model in the same way we do the data because the 
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Figure 8. Potential metal-poor moving groups in the fields at I = +10° , —20° that are located within Rq < 4.0 kpc. 



model has no metallicity or age information. Furthermore the com- 
ponents are overlapping; stars that we identify with the boxy/peanut 
bulge (components A and B) do extend in the tail down to [Fe/H] ~ 
-1.0 (see Figure[T). The single disk model is, however, unlikely to 
include the thick disk stars that we associate with component C in 
the same proportion as in the Galaxy. In Figure[TT] the observed and 
model velocities and dispersions are a fairly good match, although 
the rotation of the model is faster than the rotation of the stars and 
the dispersion of the model in the central regions at 6 = —5° is 
about 15% larger than the data. At higher latitudes, the dispersion 
flattens in both the model and data although the dispersion of the 
model is about 10% too low at |/| > 10°. The model has again been 
scaled in amplitude to best match the dispersion at 6 = —7.5° out 
to 1^1 = 20°. Note that the velocity scaling factors for the model 
used in Figures [Tol and [TTI are different, and were adjusted to best 
fit each data set for stars [Fe/H] > -0.5 and stars [Fe/H] > -1.0, 
respectively. A larger scaling factor for the velocity code units is 
used to best match the dispersion for stars with [Fe/H] > -1.0. The 
bar is rotated at about 25° to best fit the dispersion profile for both 
Figures. 

As already mentioned, our model has a single disk and no star 
formation and this could explain why its final bulge profile does 
not reflect all the populations we see in the central region of the 
Galaxy. Models that have two initial disks may offer more insight. 
If component C, the thick disk, is a separate structure, it may influ- 
ence the dynamical instability mapping of stars from the original 
disk into the bulge without itself being mapped into the X-shaped 
structure. Furthermore, given that we find two distinct components 
A and B for the bulge, and associate A with the colder thin disk in 
the inner region, it may not be surprising that we can not find an 
exact match between model and data. The origin of the duality we 
see in the bulge, with a kinematically colder metal-rich component 
and kinematically hotter more metal-poor component extending to 
higher latitudes, is not clear. It may be due to the different redistri- 
bution of colder higher metallicity disk stars and lower metallicity 
disk stars with higher orbital energies (as seen in the solar neigh- 
bourhood) to larger heights from the plane during the instability 
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Figure 9. Surface density projections of our N-body model scaled to dimen- 
sions of the Milky Way. The Sun would be located at x=0, y = -8 kpc and 
the bar is at an angle of about 20° to the Sun-centre line. The two groups of 
individual stars plotted correspond to the two groups shown in Figure[8]at 
positive (near side of the bar) and negative (far side of the bar) longitudes, 
respectively. 



event in a single disk scenario. Alternatively it may be a conse- 
quence of the instability of two disks or else a second star formation 
event during the instability process close to the plane. Addition- 
ally, although the model was scaled to the light profile seen in the 
COBE image, the relative scales and densities of the components 
of thin disk, bulge and thick disk in the bulge regions will affect 
the dispersion profile. It seems unlikely that a simple single-disk 
model in which all these stars have formed before the instability 
event will match exactly over all {l,b). Adjusting the angle of the 
bulge with respect to the Sun within a wider estimated range of 
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Figure 10. The rotation and dispersion of our model (lines) plotted with 
our data (points) for stars with [Fe/H] > -0.5. The red diamonds are for 
b = —5°, the yellow stars are for b = —7.5° and the blue rectangles for 
b = —10°. The data errors in the vertical direction ai'e smaller than the 
symbols and the error bars in the horizontal direction represent our two- 
degree fields). 




5 10 15 20 



Galactic Longitude (deg) 

Figure 11. The rotation and dispersion of our model (lines) plotted with 
our data (points) for stars with [Fe/H] > -1.0. The red diamonds are for 
b = —5°, the yellow stars are for b = —7.5° and the blue rectangles for 
b = —10°. (The data errors in the vertical direction are smaller than the 
symbols and the en'or bars in the horizontal direction represent our two- 
degree fields.) Note the scaling from code units to km s~^ is different to 
that of Figure [Tol 

15° to 30 ° ('Stanek et al.ll 19971 : iBinnev. Gerhard. & Spergelll 19971 
lBissantz"& Gerhard 2002) did not significantly improve the match. 
Further work will be done using N-body models and adjusting the 
input parameters to fit the observational data to better constrain the 
models. 



8 DISCUSSION 

In Section 7 we briefly described the comparison of our data to 
a sequence of models with different halo core radii: the best fits 
were found for models with large core radii. These will have max- 
imum disk rotation curves, while those with the small halo core 



radii will be initially clearl y sub-maximum and will evolve towards 
maximum with time (see lAthanassoulall2012l for a review). Un- 
fortunately, this finding does not allow us to set constraints on the 
dark matter distribution in the central parts of our Galaxy, because 
our models do not include a gaseous compo nent. As discussed in 
lAthanassoula. Machado & Rodionovl ^20121) . adding such a com- 
ponent will lead to a less strong bar and thus to a different quality 
fit. In other words, the properties of the bar do not depend only on 
the halo mass distribution, but also on the gas, its evolution and its 
properties. We thus cannot set any set any strong constrains on the 
former without having a reasonably solid knowledge of the latter, 
including its evolution in the past. 

We see from the ARGOS data that the boxy/peanut bulge pop- 
ulation B is mostly the fastest rotating of our bulge components (Ta- 
ble[2]l. We associate the metal-rich component A chemically with 
the younger more metal-rich stars of the inner thin disk, and its ro- 
tation and dispersion is typically ^ 80% that of component B. This 
difference in rotation between A and B is a potentially important 
signature that we can look for in chemodynamical or age-tagged 
N-body models of the instability event. In the solar neighbour- 
hood, the velocity dispersion of the thin disk is relatively small. 
However, the thin disk d ispersion decreases with Galactic radius 
jLewis & Freemanlll989l) and all components, including the most 
metal-rich boxy/peanut-bulge population A that we associate with 
the thin disk, are kinematically hot in the inner regions. 

We now discuss the distribution of velocity dispersion for the 
components of the bulge. At 6 = —10°, stars with [Fe/H] > — 1.0 
show a fairly flat distribution of velocity dispersion with longitude 
for \l\ < 20°. The velocity dispersion of components A and B 
shows structure at & = —10° at longitudes of —5°/ + 10° degrees 
(see Figure |6) which correspond approximately to equal distances 
of about 2 kpc along the inclined bulge. The velocity dispersion is 
higher in these fields, particularly for component A. This structure, 
which is not seen in the more metal-poor component C, must cor- 
respond to some feature of the orbit structure in the bar/bulge, but 
it is not apparent in our N-body models and we cannot at this stage 
identify its nature. The Sagittarius dwarf stream is present in this 
direction, but it is unlikely that we are including members of this 
stream in our sample, because we select stars within y = ±3.5 kpc 
of the centre of the Galaxy. However, the effect of any contami- 
nation from the Sagittarius dwarf stream would be to cause a fake 
high dispersion in some fields, particularly at \l\ — 5° , |6| = 10°. 
Note all the dispersions are calculated from the standard deviation 
of all stars in the sample, without a-clipping. 

For stars with [Fe/H] > —0.5 at 6 = —5°, the velocity dis- 
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persion decreases almost line arly with longitude at abo ut —12 km 
s~^ kpc^^. For comparison. iLewis & FreemanI ( Il989h measured 
the radial component of the velocity dispersion for the stars of the 
old disk of the Milky Way, from Rq = 18 kpc to about 3 kpc, 
and also the azimuthal component of the dispersion over a more 
restricted range in Ro ■ In the inner region, they find that the disper- 
sion component decreases with increasing radius at about —12 
km s^^ kpc^^, very similar to the gradient in dispersion observed 
for our c omponents A and B. At t he value of Re, corresponding to 
I = 20°. lLewis & FreemanI ( fl98^ find that is about 75 km s , 
a little higher than observed for components A and B (see Figure 
5). We note that the lLewis & FreemanI ( Il989l) field is at a latitude 
of fe = —4°, closer to the plane than our h — —5° fields. 

In our b = —5° fields, our obsei^ved stars are predominantly 
boxy/peanut-bulge members (A and B) with a smaller fraction of 
thick disk stars (C). The velocity dispersion gradient measured 
from our survey for these c omponents A and B is in very good 
agreement with the results of iLewis & FreemanI (1 1989!) for the old 
disk. The similar gradients for A and B suggest these populations 
are dynamically associated, and this would be consistent with bulge 
formation out of the disk. We see that component A has slower 
rotation and smaller velocity dispersion than component B. From 
Jeans' equation, that means that the scale length of A is shorter 
than for component B. We also argued earlier that A is more con- 
centrated to the Galactic plane, so it has a smaller scale height than 

B. It seems that A is a more compact version of B, in both verti- 
cal height and radius. For component C, the observed decrease in 
velocity dispersion with longitude is much slower (about —7 km 
s-i kpc~^) at both b = -5° and b = -10°, from I = 0° to 
20°, corresponding approximately to Rq = to 3.0 kpc. We asso- 
ciate component C with the thick disk. The different structure and 
kinematical properties of component C suggest that its dynamical 
history has been different from that of boxy/peanut bulge itself. 

There is another scenario for interpreting components A,B and 

C. Using an N-body model, it is possible to show that an underly- 
ing classical bulge can absorb angular momentum from a form- 

J " k 

ing bar dAthana ssoula 2003]). This will transform an initially small 

isotropic non-rotatin g classical bulge into a triaxial and cylindri - 
cally rotating object dSaha. Martinez- Valpuesta. & Gerhardll2012h . 
This object is then rotating along with the bulge formed from the 
disk via the dynamical instability. This transformed structure could 
be seen as one of our metallicity components. Could this structure 
be our component B or C ? This depends on the detailed predictions 
of the model: is the transfoimed classical bulge predicted to show 
the peanut structure associated with a split in the red clump stars, 
and how does its expected velocity dispersion compare to that of 
the instability-generated bulge? Given the characteristic dispersion 
seen in B, however, it seems more likely that component C would 
be a candidate for this scenario. We cannot discuss this scenario 
further wit hout more predictions from the simula t ions. An expan- 
sion of the lSaha. Martinez- Valpuesta. & GerhardI 1 I20T2I) results to 
evaluate the kinematics and spatial distribution of the inner bulge 
and the transformed classical component would enable further in- 
vestigation of this formation scenario in the light of our data. 

The rotation curves for the metal-poor stars with [Fe/H] ^ 
— 1.0 are typical of a slowly rotating population and are probably 
part of the halo and metal-weak thick disk populations. Their ve- 
locity dispersions are typically 100 — 150 km s~^, consistent with 
identification as halo stars. 

We find a kinematically distinct population with [Fe/H] < - 
1.0 which is not aligned with the main body of rotation of the 
stars in the bulge region. These stars are located in our fields 



at {l,b) = (10°, -7.5°) and (/,6) = (10°, -10°) and have 
mean metallicities ~ —1.7. These stars may belong to a stellar 
stream associated with an accretion event of a small system by the 
Galaxy. There is another possibility for the origin of the very metal- 
poor stars observed in the inner Galaxy in our survey. Accord- 
ing to ACDM simulations, the bulge region between longitudes of 
10° to 20° is expected to contain the oldest stars in the Galaxy, 
likely only a few generations older than the first stars (iTumlinsonl 
I2OI0I) . These stars are not necessarily the most metal-poor stars 
in the Galaxy but they are formed at early times in the high den- 
sity pre-galactic fragments that are rapidly chemically enriched 
due to a fast star formation rate and are subsequently accreted by 
the assembling Galaxy. Stars with [Fe/H] < —1.0 also include 
the potential first generation star candidates now located in the 
bulge (W hite & Springel 2000; Diemand, Madau, & Moore 2003; 
lTumlinsorJl2010l). These stars are expected to show distinct chemi- 
cal markers ( Kobavashi, Karakas. & Umedall201ll) and be on more 
tightl y bound orbits than ordinary stars of the halo (iTumlinsonl 
These stars may now be associated in metallicity but are 
unlikely to be associated in kinematics due to the short mixing 
timescales in the central bulge. More detailed chemical abundance 
analysis will reveal more clues as to the origin of the group of seven 
metal-poor stars with similar radial velocities and metallicites in 
our fields at {l,b) = (10°, -7.5°) and {l,b) = (10°, -10°). It 
seems likely however, that these stars may belong to more conven- 
tional stellar streams associated with small systems accreted by the 
Galaxy (see Figure [§). A more detailed chemical analysis of all of 
our stars with [Fe/H] < -1.0 will similarly test if these stars are 
all consistent with halo stars in the inner region, with comparable 
properties to the halo stars observed near the Sun, or if chemical 
markers indicate they belong to the first generation of stars pre- 
dicted by ACDM simulations. 



Comparison of the stellar kinematics as a function of the 
metallicity across longitude at b = —5° and b = —10° shows 
that stars with [Fe/H] > —1.0 have similar kinematical properties, 
which are dissimilar from th e rota tion of stars wit h [Fe/H] < —1.0. 
iHarding & Morrison! ( Il993h and iMinnitil d 19961) reported such a 
kinematic break at [Fe/H] « —1.0 from surveys of K giants in 
the bulge. They also found an increasing dispersion with decreas- 
ing [Fe/H] in their bulge fields. We find this trend of decreasing 
dispersion with metallicity from the bulge out into the disk, for 
stars with [Fe/H] > —1.0. We also see that the dispersion shows 
a discontinuity at [Fe/H] = —1.0, where both the rotation and dis- 
persion become a changing function of longitude and latitude. We 
conclude that stars with [Fe/H] > —1.0 share a common history, 
showing similar overall kinematics and a smooth transition out into 
the disk. 



Although stars with [Fe/H] > —1.0 are rotating cylindrically 
and are kinematically distinct from stars with [Fe/H] < —1.0, the 
stars with [Fe/H] > —0.5 have a clear kinematical identity. The 
stars in our components A and B show the same dispersion trends 
across (/, 6), and these are not seen for the stars with —0.5 > [Fe/H] 

> —1.0 in component C (see Figure |6j. The stars with [Fe/H] 

> —0.5 (components A and B) appear to have experienced similar 
dynamical processes which were not shared by the stars in compo- 
nent C. We will be building a model of the bu lge/bar and installin g 
it into the Galactic modelling tool Galaxia jSharma et al.ll20ll]) . 
Then we can subtract the model from the Milky Way data and look 
for substructure. 
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9 CONCLUSION 

Simulations of bar formation predict tliat tlie Galactic bulge is just 
the bar seen edge-on and not far from end-on, puffed up into a boxy 
morphology in the inner regions. We find cylindrical rotation for 
stars with [Fe/H] > —1.0 which transitions smoothly out into the 
disk. Across our three latitude zones, this rotation compares well 
with that for our N-body model of a boxy bulge formed via dynam- 
ical disk instabilities. That our N-body model does not replicate the 
velocity dispersion measured across all latitudes may reflect that 
the initial inner Galaxy probably comprised multiple populations 
with distinct kinematics and density gradients, while our simple 
model evolved from a single thin disk. 

Based on the relationship between kinematics and [Fe/H], we 
have argued that stars with [Fe/H] > — 1.0 share a common origin, 
distinct from metal-poor stars with [Fe/H] < —1.0. The character- 
istic dispersion profiles seen for stars with [Fe/H] > —0.5 suggest 
that these stars have shared a common dynamical history. The dis- 
tinct dispersion profiles for stars with [Fe/H] > —0.5 are compat- 
ible with the split stellar magnitude distribution of the red clump 
stars in the inner region which is only seen for stars with [Fe/H] 
> —0.5. The metal-poor stars with [Fe/H] < —1.0 show a more 
spheroidal rotation profile although this is not necessarily evidence 
of an underlying merger-generated bulge component. We propose 
that these stars are probably part of the halo and the metal-weak 
thick disk. We find kinematically distinct metal-poor groups of stars 
in our sample. These stars may belong to an independent counter- 
rotating stream or moving group in the inner Galaxy. 

We are presently working on more complete simulations, in- 
corporating chemical evolution and extra components, such as 
a gaseous disk, a classical bulge and a thick disk. Our obser- 
vation that the most metal-rich component A is rotating more 
slowly than component B will provide an interesting constraint 
on these more complete evolutionary chemodynamic models. One 
of the challenges for future modelling is to find a way to esti- 
mate or constrain the fraction of an underlying classical bulge 
co mponent in the Milky Way. We have s e en fro m the simulations 
of ISaha, Martinez-Valpuesta, & GerhardI (|2012|) that this can no 
longer be done reliably from the stellar kinematics alone. 
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